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Abstract 

In many-body perturbation theory (MBPT) the self-energy S = iGWT plays the key role since 
it contains all the many body effects of the system. The exact self-energy is not known; as first 
approximation one can set the vertex function T to unity which leads to the GW approximation. 
The latter properly describes the high-density regime, where screening is important; in the low- 
density regime, instead, other approximations are proposed, such as the T matrix, which describes 
multiple scattering between two particles. 

Here we combine the two approaches. Starting from the fundamental equations of MBPT we 
show how one can derive the T-matrix approximation to the self-energy in a common framework 
with GW . This allows us to elucidate several aspects of this formulation, including the origin of, 
and link between, the electron-hole and the particle-particle T matrix, the derivation of a screened 
T matrix, and the conversion of the T matrix into a vertex correction. The exactly solvable 
Hubbard molecule is used for illustration. 
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I. INTRODUCTION 



The GW approximation (GWA) to the electron self-energy S [1] is nowadays the method 
of choice for band-structure [2, 3] and photoemission calculations (see e.g. Ref. [4-7]). In 
the general expression S = iGWT it approximates the vertex function F 1, and keeps 
only the 1-particle Green's function G and the screened Coulomb interaction W. In the 
high-density regime, where screening is important, the GWA works reasonably well; in the 
low-density limit, instead, where the quantum nature of the system dominates, GW shows 
some failures [8-13]. One hence needs to go beyond GW. Iterating Hedin's equations further 
seems the obvious thing to do, but this is technically difficult and there is no guarantee that 
results will quickly improve. To overcome the problems of finite-order corrections to the 
self-energy one might then use the strategy of creating an infinite number of diagrams via 
a Dyson equation for the Green's function. This is what is done in GW, where the Green's 
function is determined to infinite order in W using a Dyson equation with a kernel (the 
self-energy) that is only linear in W. A similar strategy can be used for the self-energy itself 
by introducing the scattering T matrix [14-17] that yields the self-energy as S = GT. In the 
so-called Bethe-Goldstone approximation, originally introduced in the nuclear many-body 
problem [18] for the two-particle Green's function, T describes multiple scattering between 
two particles (two electrons or two holes) or an electron and a hole. The approximation is 
justified in the limit of low density of electrons or holes, i.e. close to completely filled or 
completely empty bands. 

T-matrix approaches have been extensively used in the context of Hubbard models [19-23] 
and the results have confirmed that this approximation is very good at low-electron density, 
precisely where GW fails, but is not superiour to the GWA at half filling where the correlation 
gap is not well reproduced [22]. However in general none of the three possibilities (GWA, 
particle-particle or electron-hole T matrix) will give an exhaustive description. This refiects 
the dilemma of how to decide which two-particle correlation to privilege in the description of 
a (at least) three-particle problem. It suggests to work with combinations. There are several 
ways to combine different correlation channels, and care must be taken to prevent double 
counting of low-order terms [24, 25]. The "fluctuating exchange" (FLEX) approximation 
by Bickers at al. [26, 27] starts from the second Born self-energy, i. e. the exact self- 
energy to second order in the Coulomb potential, and then sums all contributions, starting 
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from the third order, in each channel separately (i.e. there are no mixed diagrams). It 
contains hence the GWA and its exchange counterpart to all orders, plus, starting from third 
order, all T-matrix particle-particle diagrams, and all T-matrix electron-hole diagrams. In 
particular for the Hubbard model with its local interaction, where the exchange contribution 
simply leads to a spin- dependence of the interaction, calculations remain quite efficient. The 
Bethe-Goldstone T-matrix and GW approximations can be regarded as an approximation to 
FLEX where only one channel is taken into account. Beyond the summation of independent 
channels two-step FLEX approaches have been proposed, where one channel enters the 
calculation of a second channel through an effective screening [28-30]. One can also couple 
particle-particle and electron-hole channels on an equal footing, which then leads to the 
rather involved parquet theory [24, 31]. In the same spirit in Ref [32] a variational functional 
of the Green's function and the two-particle scattering vertex (which is a T matrix as defined 
in [33]) has been proposed. Its systematic construction yields the particle-particle T-matrix 
approach as simple approximation, and adds electron-hole diagrams at a higher level of 
perturbation theory. 

In this work we propose an alternative way to derive the coupling of GW and T-matrix 
channels on an equal footing starting from exact many-body equations. Attempts to go 
beyond GW by summing the GW self-energy to screened versions of the T-matrix self-energy 
are found in literature [12, 34-39]. However only putting the different approximations on the 
same footing one can get unambiguous corrections to GW from the T matrix. These screened 
T-matrices indeed require appropriate double counting corrections to keep the second-order 
terms exact and to avoid negative spectral functions. Therefore in the following we present 
a unified framework that links GW, GWT and T matrix. This allows us to address several 
questions, in particular. What is the origin of, and link between, the particle-particle and 
the electron-hole contributions to the T matrix? How do we get a screened version of the T 
matrix? How do we translate the physical content of the T matrix into a vertex correction? 
These questions will be answered in Sec. II. In Sec. HI we will then apply the T matrix 
to the Hubbard molecule at 1/4 and 1/2 filling. This system allows us to compare the T 
matrix, GW, and the exact results, hence to illustrate the performances of the different 
approximations. Conclusions are given in Sec. IV. 
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II. A UNIFIED FRAMEWORK 



In order to use a common language for GWT and T matrix, we start from the following 
exact expression for the self-energy: 

S(ll') = -zt;,(l+2)G2(12; 32+)G'-^(31'), (1) 

where (1) = (ri,cri,ti), (l"*") = (ri,(Ti,t^) with tf = ti + 6 {6 ^ 0+) describe space, 
spin and time coordinates, and integration over indices not present on the left is im- 
plicit throughout the paper. By adding a perturbing potential Uext and using the relation 
^2(12; 32+; [Uext])\u^^,=o = G(13)G'(22+) - _ [40], (1) can be written as 



where = -Cfg^G is used. Here vh{1) = -iv{l+2)G{22+) is the Hartree potential 
and the second term on the right-hand side defines the exchange-correlation contribution 
to the self-energy, S^-c. With the help of the Dyson equation for G, Eq. (2) can be further 
rearranged as 

S(ll') = VH{l)S{ir) + S,(ll') + zt;,(l+2)G(13)H(35; 1'4)L(42; 52+), (3) 

with S^(ll') = iv{l+V)G{ll'), S(35; 1'4) = the effective interaction, and L(42; 52) = 

„T /L the time-ordered "response" of the system to an external perturbation Uext- 

This way to write the self-energy directly displays the physics behind it, i.e. the description 
of a particle interacting with the system: the particle can scatter against the density of the 
system (Hartree term), it can exchange with another particle of the system (exchange term), 
it can do something to the system (last term), i.e. it can have an effective interaction with 
the system (H), the system responds (L), and the particle feels this response through the 
Coulomb interaction (vc). 

There are two essential ingredients in Eq. (3): the effective interaction 2(35; 1'4), and the 
response of the system L(42; 52). Combining approximations to H and to L, various approx- 
imations to the self-energy can be created. In situations where the screening is important 
one should make an effort to obtain a good L, whereas in situations where the quantum 
nature of the interaction is important [41] one would concentrate on S, although L and S 
are of course in principle linked through the Bethe-Salpeter equation [33] and one might 
wish to keep them approximately consistent. 



A. How to get GWl 



Neglecting the variation of Y^xc in i.e. keeping only the classical interaction Wc, one ob- 
tains S^c(ll') = S^ + zfc(12)G'(ll')^^c(l'4)x(42), with x(42) = -zL(42;42) the time-ordered 
response function. Hence one gets a screening contribution with respect to T^x'- this is the 
GW form, with W = Vc + VcX^c- At this stage it has not been specified yet how to calculate 
the screening: different approximations to the screening will give the various GW flavours 
(e.g. GW'^^^ and beyond [42]). If one keeps an approximate Tj^c in S one goes beyond GW 
and includes vertex corrections. For example, approximating Y^xc by the exchange-correlation 
potential of DFT, Vxc, one gets = Yx+iv,{l2)G{lV) [v,{VA) + /.e(l'4)] x(42) , where 

fxc = ^1^; this leads to S^-c = iGWT with an approximate vertex function F = 1 + fxcP-, 
where P = iGGT is the irreducible polarizability and we used x = -P + PvcX [43] . 

B. How to get the T matrix? 

One could also use the rough approximation 1/(42; 52) 
G(42)G'(25) but concentrate on a clever approximation for 
self-energy (3) as 

S(110 ^ vnimil') + S,.(ll') + zt;,(12)G'(13) 

1. An effective 4-point interaction O 

It still remains to find an appropriate approximation for the functional derivative on the 
right-hand side of Eq. (4). One can introduce an effective 4-point interaction O such that, 
similar to GW, 

E(ll') = G'(42)0(12;1'4) (5) 

Note that Eq. (5) is closely related to the expression of the self-energy within the T-matrix 
approximation as given e.g. by Kadanoff and Baym (see Eq. (56) in Ref. [17]), which 
is the goal of this derivation. However at this stage O is not yet the T matrix. Since 
G(42)0(12; 1'4) cannot be inverted to find O, several choices of O make the correct S [44]. 
First note that in (4) there are direct and exchange terms. Therefore it is convenient to 
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1. This modifies the exact 



<5S(31') 
6G{A5) 



G(42)G(25) 



(4) 



divide the self-energy as S = Si + S2 and, consequently, O = Oi + O2 with 

'5Si(31') 



Oi(12; 1'4) = -zt;,(12)(5(ll')5(42) + zi;e(12)G(13) 
02(12; 1'4) = it;e(12)5(2l')5(41) + Wc{l2)G{U) 



50(45) 
5X2(31' 



■G(25) 



G(25) 



(6) 
(7) 



5G(45) 

This decomposition of O allows us to find two interaction channels as in the T-matrix self- 
energy with a direct term (here given by Oi) and an exchange term (here given by O2) (see 
e.g. Eq. (13.23) of Ref. [45]). 

The solutions (6) and (7) are not unique; one could equally have written 



S(ll') = ^(25)0(15; 1'2) := ^(25) [Oi(15; l'2) + ©2(15; l'2)] := Si(ll') + ^2(11'), (8) 



with 



Oi(15; 1'2) = -iv^{l'2)6{ll')5{h2) + zt;,(12)G(13) 
02(15; 1'2) = iwe(l'2)5(12)(5(51') + iWe(12)G(13) 



5Si(31') 



_ 5G(45) 
5^2(31' 



G(42) 



G(42) 



(9) 
(10) 



. <^G'(45) 

These two decompositions of the self-energy are equivalent, i.e. they give the same self- 
energy if the exact S is used. 



2. A Dyson equation for O: the particle-particle and electron-hole T matrix 

At this level we do not have yet a closed expression for O, but Eq. (5) suggests an 
approximation to the functional derivative, in analogy with what one usually does in the 
framework of Bethe-Salpeter calculations based on GW [46-48] [49] : 

§11^ 0.(35; IM). (11) 

Note that, with (4) and (11), O is only an approximation to the total interaction S. Note 
also that approximation (11), together with Eq. (5) and reference [44], is along the same 
line of the approximations based on the Ward identity used, e.g., in Refs [50-53]. 

This approximation used in (6) and (7) allows one to determine O from an integral 
equation 

Of (12; 1'4) = -iVc{l2)6{ll')5{A2) + it;,(12)0(13)0(25)Of (35; l'4), (12) 
Of (12; 1'4) = iv,{l2)5{2l')5{4.l) + z?;,(12)G(13)G(25)Of (35; 1'4). (13) 



Here the subscript pp indicates the particle-particle nature of the interaction O that comes 
from the time ordering of the kernel GG (see App. A). Note that the zeroth order term in 
Eq. (13) has simply a minus sign and the indices 4 and 1', which, moreover, are external 
indices in Of^, exchanged with respect to Eq. (12) for therefore one can relate to 
Of as 

Of (12;1'4) = -Of(12;4l'). (14) 

= Of + Of with Of and Of given by (12)-(13) can be identified with the particle- 
particle T matrix [14, 45]. Using relation (14) one can verify that the pp T matrix meets 
the cross relation 0^'^(12; 1'4) = —0^^(12; 41') arising from the Pauli principle. Moreover, 
one can also verify that 0^^(12; 1'4) = 0^^(1'4; 12), i.e. the T matrix remains the same for 
a scattering process reversed. 

Equivalently, using the approximation (11) for O in Eqs (9)-(10), one also obtains Bethe- 
Salpeter-like equations 

Of (15; 1'2) = -iv,{l'2)6{ll')6{52) + it;,(12)G(13)G(42)Of (35; 1'4), (15) 
Of (15; 1'2) = tvc{l'2)6{12)6{51') + it;c(12)G(13)G(42)Of (35; 1'4), (16) 

where the subscript eh indicates the electron-hole nature of the T matrix. Note that, similar 
to the particle-particle T matrix, the zeroth order term in Of has simply a minus sign and 
the indices 1 and 5 exchanged with respect to Of. However, the fact that 1 and 5 are not 
external indices prevents from relating Of to Of, unlike in the case of the pp T matrix 
[54]. This is because the Pauli principle does not apply. The symmetry with respect to a 
reversed scattering process, instead, holds. 

Because of the approximation (11), the pp T matrix and the eh T matrix do not in 
general give the same self-energy anymore. In the first iteration of (12)-(13) and (15)-(16) 
the equality still holds, with 

EPP'(i)(ll') = S^^'«(ll') = To(ll') = VH{l)S{n') + S^.(ll') 

+ i;c(12)fc(l'5)G(ll')0(52)G(25) 

- i;e(12)^c(l'5)G(2l')0(52)G(15), (17) 

which is the second Born approximation. With the second iteration differences appear (see 
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Fig. 1 ). Indeed one obtains 

= To(ll') + Ti(ll'), (18) 
S^'''(2)(ii') = To(ll') + T2(ll') (19) 

with 

Ti(ll') = 2t;,(12)t;,(35)t;c(l'4)G(13)G(3l')G(54)G(42)G(25) 

- zt;,(12)t;,(35)i;,(l'4)G(13)G(34)G'(5lOG'(42)G'(25), (20) 
r2(ll') = zt;c(12)i;,(34)i;,(l'5)G(13)G(31')G'(54)G(42)G'(25) 

- zt;c(12)t;,(34)t;,(l'5)G(13)G(35)G'(54)G'(42)G'(2l'). (21) 

If one considers also the term 60/6G in the iteration (in (6) and (7) or (9) and (10)), then 
the equality of the self-energy via a or a eh channel is re-stablished, i.e. 

sf.f (11') = s;tf (110 = ro(ii') + Ti(ii') + T2(iiO + r3(ii') (22) 

with 

r3(ll') = it;,(12)t;,(35)t;,(l'4)G(13)G(32)G(21')G(45)G(54) 

- zt;c(12)t;e(35)i;c(l'4)G(13)G(32)G(24)G'(45)G'(5lO. (23) 

Therefore in case of the pp T matrix the two terms T2 and T3 come from the functional 
derivative 60'''^'^^^ SG , whereas in case of the eh T matrix the two terms Ti and T3 are the 
ones which are due to 60'^^'^^^ 6G. We also note that SO^^^^dG = 0, i.e. the functional 
derivative does not contribute to the first iteration, which explains why the pp and eh 
self-energy are the same after the first iteration. One can also rewrite Eq. (22) as 

Efl^^ = (SP^''(2) + E'^'i^)) - To + T3, (24) 
or alternatively, for example, as 

^% = I {^'''^'^ + ^'''^'^) + + T,) + T3. (25) 

This shows that one might try to sum the pp and eh self-energies to account for some 
terms of the functional derivatives 60 /6G, but such a sum is not obvious. For example, 
would it be worse to neglect (besides T3) the term —To or (Ti + T2)/2? We will discuss this 
issue in Sec. III. 
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Based on (3), we can now directly compare the different approximations: in GW or an 
approximate GWF, one fixes an approximation for the functional derivative of the self-energy 
to a low-order approximation, but one tries to treat L well. In the T-matrix approximation 
the latter is treated quite badly, whereas the former is kept, though approximatively, to 
infinite order in Vc thanks to the Dyson-like equation for the effective interaction. Both 
approximations allow one to describe physical processes involving three particles: the particle 
which is added to the system and the electron-hole pair that it creates. Ideally one would 
propagate the three particles together, which is not numerically affordable. Therefore one 
chooses to propagate a pair and to treat the third particle in a kind of mean-field of the 
other two. This is illustrated in Fig. (2): in GW one propagates together the electron-hole 
pair created by the additional particle, whereas in the T matrix one propagates together 
the additional electron (additional hole) and the excited electron (hole left behind) (pp T 
matrix) or the additional electron (additional hole) and the hole left behind in the electron 
excitation (excited electron) (eh T matrix). The choice is not obvious, therefore it is desirable 
to go beyond the simple scheme and to include terms coming from both the approaches. 
Before we have discussed pp and eh T matrix. Now we go towards a scheme which combines 
T matrix and GW. In order to do so we can see wether a Dyson equation for the T matrix 
can also be obtained with an improved approximation for L. 

3. Screened T matrix I 

Equations (12)-(13) and (15)-(16) define the equations for the T matrix. One can go 
beyond this approximation and include the Hartree potential or even a local part of the 
xc self-energy S^°^ in the variation 5G~^/5Ue-^f This yields L(42; 52) ^ G{A7)e~\72)G{75) 
where is a test charge-test charge screening function = 1 + vx if only the Hartree 
part is included, otherwise at least a partially test charge-test electron screening e^^ = 
X [55]. This makes the expression much more, though not fully, consistent: 
now e"^ contains a large part of the derivative of the self-energy, that is also considered in the 
effective interaction, and L contains the screening of the formerly independent propagators 
GG, that is itself based on the two-particle correlation function. Eq. (3) then becomes 

S(ll') ^vh + ^x + M'c(12)e-^(72)G(13) 



V + 



<5Ei: 



5p 



5S(310 
(5G(45) 



G(47)G(75) 



(26) 
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Note that this is the same equation as (4), but now with W = e replacing the bare 
Coulomb interaction under the integral. This leads to a screened matrix O: S(ll') = 
G(42)C)f (12; 1'4), with Of = Of^ + Ofs and 

Of^i(12; 1'4) = -iv,{12)6{ll')5{A2) + ^1^(12)0(13)0(25)0^1(35; 1'4), (27) 
0^2(12; 1'4) = = iv,{12)6{U)6{l'2) + ziy(12)0(13)0(25)Of2(35; 1'4). (28) 

In principle one could include also a nonlocal part of the xc self-energy, but the equations 
become more involved. 

As for the T matrix, also for the screened version one has 0^''^(12; 1'4) = —0^^(12; 41'). 
Similarly, one can derive the electron-hole screened T matrix, which looks like Eqs (9)-(10) 
with W replacing the bare Coulomb potential in the last term on the right-hand side. When 
= 1, and hence W = Vc, this version of the screened T matrix reduces to the one of 
Ref. [17] [56]. Note that the Hartree and exchange parts remain unscreened, which is in net 
contrast with other versions of the screened T matrix reported in literature [12, 34-39]. 

It is interesting to take the screened T-matrix equation in its first iteration: in this case 
the self-energy (particle-particle version) becomes 

Sfi«(ll') = 6{n')vH{l) + ^W{12)G{n')v,{l'4)Lo{24), (29) 
El';f\ir) = S,.(110 - W^(12)0(14)t;e(4l')0(2l')0(42), (30) 

with Lo(24) = — zO(24)0(42). The electron-hole screened T matrix produces the same 
self-energy as the particle-particle screened T matrix in its first iteration, as in case of the 
T-matrix approximation. Moreover the resulting self-energy is exact to second-order in the 
Coulomb interaction. 

In the RPA, Vr(12)wc(l'4)Lo(24) = [vc + VcLqvJ{1 - Lq)vc]LqVc = W -Vc. The sum of 
El and S2 yields hence GW, plus the last term of S2. In Fig. 3 we report the diagrammatic 
representation of this self-energy: the first two diagrams represent the Hartree and GW 
contributions, respectively, whereas the last one is a term corresponding to the second-order 
screened exchange (SOSEX). The latter contribution is becoming popular as correction to 
RPA in order to produce accurate results in the description of electronic correlation in atoms 
and solids [57]. 

We hence can conclude that GW is contained in this screened T-matrix approach, which 
moreover contains promising higher-order terms. In literature other versions of the screened 

10 



T matrix [12, 34-39] are proposed which are combined with the GW approximation to get 
the total self-energy. However, since some of the terms in the T matrix are already contained 
in the GW approximation, care must be taken to avoid double counting. In our formulation, 
instead, the screened T matrix naturally contains GW; there is hence no need to add ad hoc 
corrections. 

Because of the appearance of both Vc and W the screened T-matrix approximation (27)- 
(28) does not fulfill some symmetry conditions to be fully conserving, unlike self-consistent 
GW and (unscreened) T matrix [17]. For example, the momentum conservation law is 
violated. 



4- Screened T matrix II 

Starting from (3), one can also write 
S,,(ll') = S:.(ll') + ^fc(12)G'(13) 



svHim^i') . ^s..(3io ■ 



= T.aw + ^^e(12)G(13)^^^L(42; 52). (31) 
Using now again the approximation L(42; 52) ^ G'(47)e~^(72)G'(75) one obtains 

S.,(ll') ^ T.GW + ziy(17)G(13) ^^(45) g(47)G(75). (32) 

Note that here only the xc part of the self-energy appears. The ansatz will now 
use a screened xc matrix O: T.^^{lV) = G(42)0^p,(12; 1'4) (or equivalently S^.c(ll') = 
G(25)Og^^(15; 1'2)); the structure of the resulting equation 

0?^,(12; 1'4) = ^H^(12)5(21')5(41) + ^1^(12)G(13) [0^^,(35; l'4)G(25)] (33) 

is equal to the one determining the screened 0^\, Eq. (28), which for e^^ = 1 reduces to 
Of^. Similarly the structure of the electron-hole T-matrix Of^^(15; 1'2) is equal to the one 
determining the screened 0^25 which for = 1 reduces to Of^. Note that this screened T 
matrix yields a self-energy corresponding to the exchange term only of Refs [12, 35], where 
a screened T matrix is also used. 

Iterating this equation for O^p^ once, one gets the GWT^^^ [1, 58] approximation, namely: 

0^'j^\l2- 1'4) ^ iW{12)6{21')6{41) - iy(12)G(14)iy(4lOG'(21'), (34) 
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from which S^c(ll') = S^]^ - G {i2)W (12) G{U)W {41') G {21'). One can verify that the 



electron-hole screened T matrix produces the same self-energy as the particle-particle 
screened T matrix in its first iteration, as in case of the T-matrix approximation. It is, 
not surprisingly, also what one obtains by iterating Hedin's equations [1]: starting from 
Sj-c = iGW the equation for the irreducible vertex function reads 



To first order one gets hence 1(12; 3) = 5(12)5(13) + iW {21) G {13) G {32) and 

S^e(ll') = iG(14)r(4l'; 2)1^(21) ^ S^J^ - G{U)W{l'A)G{42)G{2l')W{21), (36) 



which is exactly the xc self-energy obtained using (34). There is instead no equivalent to 
(27) or its eh counterpart. However, in this second version of the screened T matrix, Hartree 
and exchange are not treated on the same level, since the latter is included in the definition 
of the T matrix, and the former not. This also explains why for e^^ = 1 this screened T 
matrix does not reduce to the T matrix; in fact it reduces to only a part of it, i.e., to O2. 
In Sec. Ill we will show that this unbalance prevents this version of the screened T matrix 
to give the exact result for one electron in the atomic limit, unlike the screened T matrix 



5. Vertex corrections from T matrix 

Above we have shown that the pp screened T matrix (27)- (28) (as well as the eh screened 
T matrix) contains GW plus extra terms beyond GW. We can therefore use the T matrix 
to formulate vertex corrections beyond GW. Starting from the self-energy S = GO^^ = 



G[OZ + with and given by Eqs (27)-(28), 

S(ll') ^ vh + ^^(13) [v^{13)S{l'3) + iy(12)0f (35; 1'4)G'(42)G'(25)] . (37) 
Using v^{13) = W{13) - W{12)P{25)vc{53) we can rewrite (37) as 

E(ll') ^vh + iG{13) [W{12)6{1'3)6{23) - W{12)P{25)v^{53)6{l'3) 



r(12;3) = 6{12)6{13) + iW{21)G{16)G{72)T{67;3). 



(35) 



given by Eqs (27)- (28). 



+ Vr(12)0f (35; 1'4)G'(42)G(25)] 
= Vh + iG{13)W{12)T{31';2), 



(38) 
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where 

r(3l'; 2) = 5(31')5(21') - P{25)vc{53)6{31') + Of (35; l'4)G(42)G(25). (39) 

This is an approximate vertex function that will generate the same self-energy as the pp 
screened T matrix. In an analogous way one can obtain an approximate vertex function 
that will yield the same self-energy as the eh screened T matrix: the equation will be the 
same as (39) with O^^ replacing Of on the right-hand side. It is interesting to compare this 
vertex function with the exact expression derived by Bruneval et al. [55] that reads 

r(31'; 2) = 5(31')5(21') + ^^rffV (52). (40) 

One can approximate the second term on the right-hand side of (40) as 
5S.c(31')^^52) = iM^ll^^i£(M^p(52) 



(5p(5) ^ ' 5p(5) 
(5E(31') (50(45) 



P(25)t;e(53)5(3l') 



5G(45) 

^ 0,(35; 1'4)G(42)G(25) - P(25)t;,(53)<5(31') (41) 

with Vtot = Uext + vh as the total classical potential. The last line is an approximation 
obtained with ^ ^ Os, where O, is either the pp or eh screened T matrix, and ~ GG. 
Equation (41) with (40) yields (39). It then becomes clear that the term —Pvc in Eq. 
(39), that is created by the induced Hartree potential felt by a classical particle, needs to 
be subtracted since it is already described by the GW approximation for F = 1, whereas 
the remaining part is responsible for the induced exchange-correlation potentials felt by a 
fermion. We also note that the term Pvc, which comes out in a natural way in our derivation, 
gives rise to the correlation part of the GW self-energy, = iGv^PW = iG{W — Vc). It 
equals the correction done a posteriori in the other existent formulations of the screened T 
matrix to avoid double counting [12, 35, 59]. 



III. APPLICATION TO THE HUBBARD MOLECULE 

In the following the performance of the T matrix as well as the screened T matrix I, 
and hence of the approximate vertex function (39), is illustrated using the exactly solvable 
Hubbard molecule at 1/4 and 1/2 filling [60, 61] discussed e.g. in Ref. [13]. The Hubbard 
model is traditionally used to model strongly correlated systems; these are precisely the 
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systems for which GW shows failures. This, together with an exact analytical solution at 
hand, represents a powerful tool to test the improvements over GW due to the inclusion of 
pp and eh correlation for low-density systems. 

The Hamiltonian of the Hubbard molecule reads 

ij = l,2 a i=l,2 aa' a,i=l,2 

Here njo- = cl^Ci„, and Qo- are the creation and annihilation operators for an electron 
at site i with spin a, U is the on-site (spin-independent) interaction, —t is the hopping 
kinetic energy, and eo is the orbital energy. The Hamiltonian further contains a potential 
Vo that can be chosen to fix the zero energy scale. The physics of the Hubbard model 
arises from the competition between the hopping term, which prefers to delocalize electrons, 
and the on-site interaction, which favours localization. The ratio U/t is a measure for the 
relative contribution of both terms and is the intrinsic, dimensionless coupling constant of 
the Hubbard model, which will be used in the following. 

Projected onto the (orthonormal) site-basis of the Hubbard model the T matrix (12)-(13) 
and (15)-(16) becomes 

07^'{^) = -^SaS,, [0-7(0;) - 5.,.,0--(a;)] , (43) 

with O^'^^iuj) = [1 + ULl''"'{uj)]'^U, from which 

Kjioo) = -^J '^^G%{u)0:r{uj±u). (44) 

Here a indicates a spin opposite to a, the sign '+' refers to the particle-particle contribution 
for which Lq^J^'^^{u}) = —i J ^G^j{uj')G^j{uj — uj'), and the sign ' refers to the electron- 
hole contribution for which ^^^."^•"'^(w) = -i J '^G1]{u')G1]{u' - u). More details on the 
spin-structure of the T matrix can be found in App. B. 

The equations for the screened T matrix are more involved: because the screened Coulomb 
interaction W is nonlocal in space and frequency-dependent (see Ref.[13] and C7), where 
a RPA W is given), one has to solve a 4-point equation similar to the dynamical Bethe- 
Salpeter equation for electron- hole excitations [47, 62, 63]. This is beyond the scope of this 
paper, and of the majority of the applications. Moreover, we observe that when there is no 
screening in the system (which is the case in the atomic limit t — ?■ for the model used here) 
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the screened T matrix reduces to the T matrix, and one hence retrieves the same results 
as those obtained with the latter. For a finite t we assume that for the model used here 
the major contribution to the T matrix arises from the on-site screened interaction. This 
is dominated by the bare interaction U, which justifies to take the screened interaction in 
its static {u = 0) limit. In this case the structure of the screened T matrix is the same as 
for the T matrix with the onsite screened Coulomb interaction W = U ~ (i±^ivp2)c^^ with 
^2 _ -\-2Ut{l + 6N=2) and the total number of electrons in the system (see Ref.[13] and 
C7)), replacing U. We notice that with this approximation, in particular assuming a static 
W, the screened T matrix will not reduce anymore to the T matrix in the atomic limit. This 
is because in the atomic limit the frequency- dependent part of the screened interaction, i.e. 
■^i^^j^rfl^y^! vanishes if w 7^ 0, but will reduce to — y if w = 0. 

Another simplification that will be adopted in the following is the neglect of self- 
consistency. In principle the T-matrix approximation requires self-consistency in order 
to be conserving. The various levels of self-consistency in the T matrix have been ad- 
dressed by several authors in particular in model systems (see e.g. [21-23, 64-66]). When 
the quantity of interest is the spectral function, as in our case, self-consistency was found 
to deteriorate spectra [21-23, 67-69]. In practice, as in the case of the GWA, rather a 
nonself-consistent "best G" strategy is adopted. 



A. Hubbard molecule 1/4 filling 

We first consider the Hubbard model with only one electron in the ground state; the 
ground state is |\E'o) = l/\/2{\ '\ 0) + |0 f))) where the electron has spin up (equivalently 
the spin-down situation could be chosen). 



1. T matrix 

Using the noninteracting Green's function for 1/4 filling of Ref. [13] [Eqs (16) and (18)] 
and Eq. (44) the particle-particle contribution to the T matrix yields the exact self-energy, 
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which reads 



= 



f/2 



1 



-1) 



(45) 
(46) 



w — Cq — -J — 3t + i?7 u — eQ — ^ — t + ir] 

This self-energy gives rise to the exact one-particle Green's function and spectral function. 

The eh T matrix also yields the exact result for the spin-up self-energy, but a spin-down 
self-energy with poles that are shifted by U from the exact ones: 





u 



2 + 



f/2 



+ 



CO - eQ + ^ - 3t + irj u-eQ-t + ^+ir] 



(47) 
(48) 



The poles of this spin-down self-energy give addition and removal energies in overall bad 
agreement with the exact result, as shown in Fig. 4, where exact, T-matrix, and GW renor- 
malized addition energies u"'/t for the spin-down channel are reported versus U/t. Par- 
ticularly interesting is the spectral function at U/t — ?■ oo: as discussed in Ref. [13], when 
t — (atomic limit) the electron is localized either on one site or the other with the same 
probability. Therefore two electron peaks, with the same weight 1/2, appear in the spectral 
function (see Fig. 4), one corresponding to the addition of the second electron on the empty 
site (peak at eg) and the other corresponding to the addition on the filled site (peak at 
^0 + U). The GWA produces only one peak at + U/2 with spectral weight 1. This is 
due to the interpretation of the charge density as an average charge distribution rather than 
a probability. The T matrix, instead, "sees" where the electron is, although only the pp 
T matrix "sees" well. Indeed, the eh T matrix yields two peaks with the correct spectral 
weight, but at the wrong position, namely eo + U\^/2 and eQ — U -\/2/2. Therefore, it is clear 
that in the case of the Hubbard molecule with one electron the particle-particle contribution 
to the T matrix describes the essential physics. 

Note that the eh T matrix at first iteration performs as the pp T matrix, since S'^'^'^^^ = 
j]pp,{i)^ as we have already shown. In the atomic limit the T matrix at the first iteration 
shows two peaks in the spin-down spectral function: one located at = eo + U{1 — v^)/4 
with spectral weight (1 — 1/^/5)/ 2 ^ 0.28 and the other one located at u = eQ + U{l + a/5)/4 
with spectral weight (1 + l/V5)/2 ^ 0.72 (see Fig. 4). It hence contains the right physics, 
although the results are still poor. 
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2. Screened T matrix 



As illustrated above, the pp T matrix yields the exact result for the Hubbard molecule at 
1/4 filling. However, if one is interested in many-electron systems, where screening becomes 
important, the screened T matrix is more appropriate. It is interesting to check how well the 
screened T matrix performs in the 1 electron limit. In order to do so we concentrate on the pp 
screened T matrix only. Using the onsite and instantaneous approximation W = U — U'^t/h'^ 
(with = At^ + 2U t) for the screened Coulomb interaction, the screened pp T matrix yields 
the self-energy 

= = (49) 



1 (-1)(*-J) 



w — Cq — -y — 3t + 2?7 u — eo — ^ — t + tf] 



(50) 



which gives rise to the spin-down addition energies reported in Fig. 5. The approximate 
screened T matrix performs in general much better than GW; in particular, in the atomic 
limit, although it does not reproduce the exact result, it produces the correct number of 
peaks in the spectral function, unlike the GW approximations which yields only one peak. 
Therefore, already in this approximate version the screened T matrix 1 contains the essential 
interaction processes. Note that the dynamically screened T matrix I reproduces exactly 
the atomic limit for the present problem. 

The screened T matrix II, instead, does not reproduce correctly the atomic limit. Indeed, 
with W = Vc the screened T matrix II reduces to the O2 component only of the T matrix and 
this is not sufficient to capture the interactions in the atomic limit. One finds the self-energy 

Sff (^) = Sj'^(u;) = (51) 
Sff ^(^) = ^t\^) = (52) 

which is the same as the one obtained within GW . 



B. Hubbard molecule 1/2 filling 



We now consider the case with two electrons in the ground state. Technical details are 
given in App. C (see also Ref. [70]). 
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1. T matrix 



In this case neither the pp nor the eh T matrix reproduce the exact result, with the 
self-energy reading as 



U . UH 



2 Ah 



1 



+ 



u! — t — h + ii] u! + t + h — iri 



(53) 



with h'^ = + 2tU, in the particle-particle T-matrix approximation, and 



(_l)(i-i) 



(54) 



oj — t — h' + ir] oj + t + h' — ir] \ 

with h'"^ = — 2tU, in the electron-hole T-matrix approximation. The pp T matrix 
performs rather well over a wide U/t range as one can see in Fig. 6, where the renormalized 
addition/removal energies u/t are plotted versus U/t. In particular the satellite energies 
(outer energies) are better described than in GW, in line with previous findings [12, 21, 22]. 
The energies calculated using the electron-hole T matrix, instead, show divergencies. In the 
— > oo limit, all approximations studied are rather poor. 

In the atomic limit there are no double occupancies, therefore the two electrons, one with 
spin up and the other with spin down, are localized one on one site and the other on the other 
site with equal probability, i.e. the ground state is the singlet l^'o) = :^(| t i) ~ I it))- The 
spectral function thus shows, for each spin, two peaks with the same spectral weight 1/2, 
one for the removal of an electron (peak at eo), and one for the addition of a second electron 
(peak at eQ + U), as shown in Fig. 6. In App. C we show that, using the noninteracting 
Green's function, GW fails also in the case of 1/2 filling, producing for each spin only one 
kind of peak at eo + f//2 with spectral weight 1/2, both for electron removal and addition. 
This can again be understood considering that GW treats the charge/spin density as a 
classical distribution, namely half electron with half spin up and half electron with half spin 
down on each atom that respond to the additional electron or hole in the atomic limit. We 
find that the particle-particle T matrix yields the same result as GW in the atomic limit, 
whereas the electron-hole T matrix shows divergencies. 

Why is the pp T matrix exact for one electron in the atomic limit, and not for two 
electrons? To derive the T matrix we used the approximation ~ GG. In the case of 
one electron this is not an approximation, but it is the exact time-ordered response, and 
therefore the (pp) T matrix yields the exact result for one electron. This is not the case for 
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two electrons for which ^ GG is a rough approximation, and one needs to include some 
screening. The screened T matrix I indeed improves over the T matrix for two electrons even 
with an approximate RPA screening, as it is shown in the next section; such approximate 
screening is instead dramatic for one electron and a more accurate screening is needed (as 
pointed out above the exact screening would yield the T matrix and hence an exact result 
for one electron). One should hence use a screened T matrix with a screened interaction 
adapted to the system. 

Interestingly the first iteration for both the particle-particle and electron-hole contribu- 
tions to the T matrix gives the exact results for all t. Indeed after the first iteration the T 
matrix reads 

[U5,,-U'Ll]^{uj)\ (55) 



with Lq being 



u — 2t + ir] u + 2t — ir] 



(_1)CW) 



1 



(56) 
(57) 



oj — 2t + irj oo -\- 2t — irj 

for the particle-particle and electron-hole contribution, respectively. The self-energy hence 
becomes 

f/ f/2 r 1 (-i)(^--') 



00 



■le/i, (71,(1) 



00 



.<Ti,(l) 



(u) = 6ij— + 



LO — 3t 



irj 



+ 



irj 



(58) 



which is the exact one. This result, however, is peculiar for the Hubbard molecule. Indeed, 
the first iteration of the T matrix corresponds to the second Born approximation, which 
has already been explored on bigger Hubbard clusters for different fillings and interactions 
[21, 22]: indeed it does not generate the exact result and in general the T matrix is the most 
accurate at low densities; only at half-filling the T matrix is not superior to second-Born and 
this is in line with our findings. Within the GW approximation, if one considers W at first 
iteration {wI^\(jo) = USij + U"^ J2a -^o resulting self-energy is not exact, but very 



close to the exact one, differing only in the prefactor of the frequency-dependent part that 
is [/^/4 instead of f/^/8. The addition and removal energies are thus improved with respect 
to GW, although the agreement with the exact result worsens with increasing U. It is worth 
noticing that if one considers also the exchange counterpart in GW, i.e. if one includes 
not only the Hartree potential but also the exchange self-energy in the self-energy variation 



19 



STj/SG in Eq. (3), then one obtains a GW-like self-energy with a modified spin-dependent 
screened interaction W^-" {u) = U 6ij + U J2r ar -^o,'jr('^) ^ ^o-,.o-J Wrj{uj) that at second order 
in U produces the exact self-energy. Indeed this is the second Born approximation. 

The exact result obtained with the T matrix at first iteration, however, deteriorates with 
the second iteration; in this case the T matrix reads 

Oip^^'\uj) = 0:r^^'\uj) + U'Y. l^{uj)V^j{uj) (59) 



and the two self-energies become 

l^) + I _ 3^ + ,^)2 - + 3t - zr/)2 " 2t 



-1) 



u — 3t + ir] u + 3t 
1 , (-1)(^ 



' ' ' ' 16 t - 3t + ir/)2 (w + 3t - 2r/)2 2t - 3t + ir/ ' to + 3t 

Combining the two interaction channels by adding Eqs (60) and (61), the second terms 
on the right-hand side of the two equations cancel each other, thus restoring the exact 
result if | (J^pP'C^) + j^^hX^)^ jg ^^ken. We have already shown in Eq. (25) that the sum 
1 (|5]PP:(2) _|_ j^eh,{2)^ takes into account some of the terms which would appear in the self- 
energy if also the functional derivative 60 /6G were considered. This might justify the exact 
result that is obtained by taking | (S^P'^^) + S^'^'^^)). Also the sum (S^P'^^) + S^'^'^^)) takes 
into account some of the terms arising from the functional derivative (see Eq. (25)); however 
it does not give the exact result. In other words it seems more important to take into 
account the term — Tq than (Ti + T2)/2 (see Eqs (24-25)) at least in the present problem. 
However with the third iteration the fourth order terms in the pp and eh T matrix self- 
energies are the same and they would not cancel each other if the sum | {YjPp^^"^) + Yj^^A"^)^ 
is taken. Instead with the forth iteration the fifth order terms in the pp and eh T matrix 
self-energies would cancel each other. In general for the present problem the pp and eh 
T matrix self-energies show, starting from second order, the same even-order terms and 
opposite odd-order terms, as one can verify Taylor expanding the frequency-dependent part 
of the pp and eh self-energies (Eqs (53) and (54)) for small U. The same holds for the 
Hubbard model at 1/4 filling. Therefore summing the two contributions will not give the 
exact result. Even adding the GW self-energy terms and its exchange counterparts, in the 
spirit of the FLEX approximation, will not produce the exact result. 

These findings show that there is no an ultimate way to combine diagrams, and this is of 
clear relevance for realistic studies where several attempts to combine pp and eh channels 
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have been done (see e.g. [28-30, 35]). A possibility is to use the screened T matrix we 
introduced in Sec. II B 3, which produces, at least for the studied problem, results overall 
better than GW and T matrix, as shown in the following. 



2. Screened T matrix 



The pp screened T matrix with the approximate onsite static Coulomb interaction W 
U - 2UH/h'^ (with /i^ = + AUt) leads to the self-energy 

1 (-1)(^^-^) 



U ^ UWt 
2 + -W 



t — h + irj uj + t + h — iri_ 

with /i^ = + 2tW . The resulting renormalized addition and removal energies u/t are 
plotted in Fig. 7 versus U /t, and compared with the exact, pp T matrix, and GW results. 
This approximate pp screened T matrix is overall superior to the GW and the T matrix 
in the selected U /t range in the left panel of Fig. 7. In the limit U/t oo the results get 
corrupted. As a consequence the screened T matrix performs as poor as GW and the T 
matrix. 



(62) 



Our findings suggest that the screened T-matrix approximation (or equivalently GWT, 
with r obtained from the screened T matrix ) is expected to describe properly also larger, 
more dense systems. This is in agreement with the idea behind the screened T matrix to 
combine T matrix and GW and to take advantage of the strength of both approaches. For 
short-range interactions, where screening is not important, the screened T matrix reduces 
to the T matrix, which is suitable for treating short-range correlation. For long-range inter- 
action, where, instead, screening is important, we find that the screened T matrix behaves 
more like GW (in its first iteration, indeed, it gives GW and SOSEX, which is actually al- 
ready used to improve GW), which is capable of taking into account long-range correlation. 
Therefore the screened T matrix is able to capture the physics of systems with effective 
short-range interactions as well as of systems with effective long-range interactions. 

For practical calculations the equations become very involved. However, the Dyson equa- 
tion for the T matrix is similar to the Bethe-Salpeter equation for the eh excitations, there- 
fore one might use similar strategies to solve it. For example one may consider the static 
approximation to the screened Coulomb interaction, as usually used for the BSE and as 
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adopted for the screened T matrix of Ref. [12]. For short-range interactions one can con- 
sider, together with the static approximation for W, also a local approximation for the 
four-point kernel GG which was proposed in Ref. [71] for a similar screened T matrix and 
validated by several applications [34, 35, 71]; this simplifies a lot the calculation of the 
screened T matrix, which becomes a two-point quantity similar to the screened Coulomb 
potential W, and opens the way for wide applications. 

IV. CONCLUSIONS 

In this paper we have given an alternative derivation of the T-matrix approximation to 
the self-energy starting from exact equations. This allowed us to: i) link the T matrix to 
Hedin's equations; ii) understand the origin of the electron-hole (eh) and the particle-particle 
(pp) T matrix; iii) derive a screened T matrix; iv) translate the physical content of the T 
matrix into a vertex correction; v) put the second-order screened exchange (SOSEX) on the 
same level as the T matrix. 

We applied the T matrix to the exactly solvable Hubbard molecule at 1/4 and 1/2 filling 
and we studied its performance with increasing ratio U/t. We found that the particle-particle 
T matrix gives the exact removal and addition energies for 1/4 filling. The electron- hole T 
matrix, instead, performs badly. 

In the case of 1/2 filling the pp performs in general better than GW, in particular in 
describing the satellite position, whereas the eh T matrix show divergencies for U/t > 2. 
In the atomic limit t — > both GW and T matrix are very poor. In their first iteration 
both particle-particle and electron- hole T matrix produce the exact result. This result gets 
corrupted with further iterations; however pp and eh self-energies have the same even-order 
terms and opposite odd-order terms (except the first order term which is the same) at least 
in the model analyzed here. This means that one gets the exact self-energy if the sum 
i ^5]p?':(2) _|_ j^eh,{2)^ ^Yie second iteration is taken, and one retains only the even-order 
terms if the eh an pp self-energies are taken at infinite order. This can be of relevance in 
realistic calculations where eh and pp channels are combined together for improving the 
results. 

We have also studied the performance of the pp screened T matrix. The screened T 
matrix I (see text), in which Hartree and exchange terms are treated on an equal footing, 
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reduces to the T matrix in the atomic hmit when a dynamically screened interaction is 
used; even with an approximate W it is better than GW at 1/4 filling, whereas at 1/2 
filling it is overall superior to both the pp T matrix and GW over a wide U/t range. This 
means that the vertex corrections derived from this version of the pp screened T matrix can 
visibly improve over GW. The screened T matrix II (see text), in which only exchange-like 
terms are included, reduces to only one part of the T matrix in the atomic limit and this 
is not sufficient to describe exactly this limit at 1/4 filling. We show that this version of 
the screened T matrix corresponds at first iteration to the GWT^^^ approximation, which 
hence is also not appropriate to treat the atomic limit. Our illustration of the different 
T-matrix approximations on the Hubbard molecule suggests that the screened T matrix I is 
a promising approximation also for realistic systems, since it combines on an equal footing 
on one hand the physics of the T-matrix approximation, which properly describes short- 
range interaction, and on the other hand the physics of GW, which is more appropriate for 
long-range interaction. 
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Appendix A: Time structure of the T matrix 

First we examine the time structure of the T matrix (12)-(13) and (15)-(16). 
1. particle-particle T matrix 

We start from Eqs (12)-(13). In the following the indices will refer to the time only. We 
can define Of (12; 1'4) := -i5(12)Of (11; 1'4) which leads to 

Of (1 - 1'; 1 - 4) = v,6{ll')6{Al+) + ^t;,G(l - 3)0f (3 - 1'; 3 - 4)G(1+ - 3). (Al) 

Here the correct order of the field operators is explicitly ensured by the infinitesimal larger 
1+ = 1 + 7/ > 1. The two Green's functions in the product have the same time order, contrary 
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to the usual electron-hole response function. This leads to a particle-particle function L^{1 — 
3) := —iG{l — 3)G(1+ — 3). We then get in frequency space 

Of (a;; u') = Vc- v,L^{uj + uj')Of{uj- J), (A2) 

with Lf{uj) = -if ^G{uj')G{u - uj^^'^K This finally imphes that Of depends only on 
the sum of frequencies: 

Of + J) =Vc- VcLl^iio + a;')Of + u). (A3) 

This is a Dyson-like equation similar to the screening equation in the GW approximation. 
The time structure of the self-energy becomes 

- 1') = G(4 - 2+)0f (12; 1'4) = -iG{A - l+)Of (1 - l'; 1 - 4) (A4) 

and hence 

Sf (w) = -i [ — G(a;')Of (w; w')e'"''' = / — G(a;')Of (w + c^Oe'"'^- (A5) 

J 2-71 J 2lT 

Again, this is very close to the structure of GW; simply, one has a "particle-particle-screened" 
interaction —Of. One can verify that the time structure of Of is the same as for Of. 

2. electron-hole T matrix 

In case of the electron-hole T matrix (15)-(16) one can do similar steps as above and 
arrives at 

Of (1 - 5; 1 - 1') = VcS{n')5{l+5) - Vc^^{l - 3)0f (3 - 5; 3 - l'), (A6) 

where we defined the electron-hole functionLQ'^(l — 3) := —iG{l — 3)G(3 — l'^). In Fourier 
space one gets 

01\lu + Lu') =v,- v,Lf{u + uj')Ot{uj + u'), (AT) 

with L^^i.uj) = -iJ ^G{u')G{u' - u)e''^''^. 
The time structure of the self-energy becomes 

Si(l - 1') = G(2 - 5+)0f (15; l'2) = -tG{l - 5+)0f (1 - 5; 1 - l') (A8) 

hence 

Ef (w) = -'I 

For Of'^ one can proceed in a similar way as above. 
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/ — G(a;')Of (u;'; w)e*"''' = / —G{u')Of{uj - ^Oe'"'"- (A9) 
J 2n J 2n 



Appendix B: Spin structure of the T matrix 



We now schematize the spin structure of the T matrix. Both for pp and eh T matrix one 
has, in the colhnear hmit, S(cr) = G{a2)0{aa2] cra2), with 

Oi(aa2; aa2) = -ivc + iVcG{o)0i{aa2] aa2)G{o2)-, (Bl) 

and 

02((T(T2; crcTs) = A<72 + iVcG{a)02{(rcr2] (ja2)G{a2). (B2) 

Note that, unhke the screened interaction W used in the GW approximation, the T matrix 
is spin-dependent. 

Appendix C: Hubbard molecule at half filling 

The starting point is the following Hubbard Hamiltonian 

H = -t ^ ^ 4cj, + ^ 5^ 5^ ct4,Q,/Q<, + eo XI + 

ij = l,2 a i=l,2 aa' a,i=l,2 

Here = c\^Cia, and Cia are the creation and annihilation operators for an electron at 
site i with spin a, U is the on-site (spin- independent) interaction, —t is the hopping kinetic 
energy, and eo is the orbital energy. The Hamiltonian further contains a potential Vq that 
can be chosen to fix the zero energy scale. The eigenstates of the system will be linear 
combinations of Slater determinants, which are denoted by the kets |1 2), with occupations 
of the sites 1, 2 given by 0, f; i, ti- We choose eo = — y and Vq = , where is the total 
number of electrons in the system, i.e. 2 in our case. This choice is particular convenient 
since the obtained Hamiltonian, 

H = -tY.Y. + ^Y.Y. c\Aa'Cia'C^a -J U,^ + U, (C2) 

j,J = l,2 a i=l,2 aa' a,i=l,2 

has high symmetry (particle-hole symmetry), as we shall see below (see also [70]). Using 
this Hamiltonian we can calculate the exact one particle Green's function. 
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1. Exact solution 



The main exact quantities we are interested in are 



2a2 



(1 



At \2 



+ 



2a? 



UJ - (c/2 -t) +iri uj - (c/2 + t) + ir] 



it \2' 



UJ + (c/2 -t) - ir] u + (c/2 + t) - irj 



(C3) 



i-j) 



u — t + irj^u + t — irj 



(C4) 



u 

~2 



6^J + 



+ 



(C5) 



u — 3t + ir] u + 3t — ir] 

with c^ = 16t^ + and = 2 (^j^z^ + 1 j • Note that the symmetry of the system is such 
that G\i = G\i = G22 = and = '-^12 = '^21 = ^21 similarly for the self-energy. 



2. GW solution 



Also here we give the main results 



4 



1 



UJ — 2t + ir] UJ + 2t — ir] 



(C6) 



W,,{uj) = US,, + (-1)(^-^') 



UH 



UJ — h + ir] UJ + h — ir] 



(C7) 



— + 



(C8) 

UJ — {t + h) + ir] ' a; + (t + /;,) — zr/J ' 
where h"^ = 4t^ + AtU. Note that only the component P'^'^ is given in (C6) (and not the full 
P°"°"'), since only this is needed to calculate W. 

The poles of the one-particle Green's function can be calculated using 



det[G-^] = det[Go^ - S] = det 



UJ — Sii t — H 



12 



t — S12 UJ — Tj 



11 



from where we get 



11 — ^12 J 



u; = -t + (S 



11 -r z^i2 



(C9) 



(CIO) 
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which is general, i.e. we can apply for any approximation to the self-energy. In the case of 
GW, we get the following poles 



Wl,2 = 



|-/i±y/(/i + f + 2t)2 + 



h 



(Cll) 



2 



t^3,4 = 



f + /i±y/(/i-f + 2t)2 + 



h 



(C12) 



2 



We note that for f/ 7^ the particle-hole symmetry is lost due to the lack of self-consistency 
using the GqWq approximation [13]. This symmetry can be enforced by absorbing the static 
part of the self-energy {U/2) into the chemical potential; this ultimately corresponds to 
dropping the terms U/2 in (C11)-(C12). 
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FIG. 1. Diagrams corresponding to the self-energy obtained with the second iteration of the 
particle-particle T matrix (S^^'^^^ and S2^'^^^)and electron-hole T matrix (S^'^'^^^ and Sj'^'^^^). 
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FIG. 2. Schematic representation of the physical contents of GW (a), pp T matrix (b), and eh T 
matrix for particles with collinear spins. 
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FIG. 3. Diagrams corresponding to the self-energy obtained with the first iteration of the screened 
T matrix. The diagrams, from left to right, represent the Hartree, GW, and second-order screened 
exchange (SOSEX) terms, respectively. Note that in the GW term we collapsed the two terms 
which for = 1, i.e. W = Vc, reduce to the second and first diagrams of '^j^l'^^ and YI^^^^ , 
respectively, of Fig. (1). 
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FIG. 4. Two-site Hubbard model at 1/4 filling: comparison between the exact spin-down renor- 
malized addition energies uj"'/t (solid lines) as function oiU/t (left panel) and Log(U/t) (right 
panel) and the results obtained from GW (dashes), particle-particle (solid lines, equal to the exact 
result), electron- hole (crosses), and 1st iteration T matrix (circles). In the atomic limit the spectral 
function, i.e. the peak positions and weights, is illustrated on the right-hand side, upon multiplying 
by t and taking the t — t- limit. 
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FIG. 5. Two-site Hubbard model at 1/4 filling: comparison between the exact spin-down renormal- 
ized addition energies to"" /t (solid lines) as function of U /t (left panel) and Log(U/t) (right panel) 
and the results obtained from GW (dashes) and particle-particle (solid lines, equal to the exact 
result) and (approximate) screened T matrix (triangles). In the atomic limit the spectral function, 
i.e. the peak positions and weights, is illustrated on the right-hand side, upon multiplying by t 
and taking the t — )• limit. 
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FIG. 6. Two-site Hubbard model at 1/2 filling: comparison between the exact renormalized 
addition/removal energies ijjjt (solid lines) as function of U jt (left panel) and Log(U/t) (right 
panel) and the results obtained from GW (dashes), particle-particle (dots), electron-hole (crosses), 
and 1st iteration T matrix (solid lines, equal to the exact result). In the atomic limit the spectral 
function, i.e. the peak positions and weights, is illustrated on the right-hand side, upon multiplying 
by t and taking the t — t- limit. 
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FIG. 7. Two-site Hubbard model at 1/2 filling: comparison between the exact renormalized 
addition/removal energies ijjjt (solid lines) as function of U jt (left panel) and Log(U/t) (right 
panel) and the results obtained from GW (dashes) and particle-particle (dots) and screened T 
matrix (triangles). In the atomic limit the spectral function, i.e. the peak positions and weights, 
is illustrated on the right-hand side, upon multiplying by t and taking the t — )• limit. 
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